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ABSTRACT 

Root locus plots in the s plane are practically limited to a span of 
two decades. In problems where some poles are very close to the origin, 
while other poles or zeros are at very large values of s, sufficient 
resolution of the locus may not be available in two decades of the s plane. 
However, root loci can be plotted in the ln s plane using as many decades 
as desired (and whichever are desired). 

Electro-Scientific Industries Corporation manufactures an analog com- 
puter, THE ESIAC, which constructs root locus plots in the ln s plane. 
This paper presents a mathematical analysis of the ln s plane, a simple 
method of plotting root loci in the ln s plane, and several practical 
examples demonstrating the technique. Construction of root loci by this 


method requires the use of one template. 
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PLOTTING IN THE LOG S PLANE 
ye Introduction 

The value of root locus plotting is limited in those applications 
where some poles and zeros are widely separated and others are close to 
the origin or close to each other. When scale is expanded to show suffi- 
cient detail close to the origin, other poles are moved off the paper. 
The designer needs to see all the poles and zeros, yet retain sufficient 
detail near the origin. 

Logarithmic scaling solves the problem. This paper examines the 
properties of the In s plane and then shows how each root locus sketching/ 
plotting technique in linear coordinates (s plane) also applies in the 
In s plane. In the s plane, all the present root locus technique is 
based upon straight lines, slopes, and lengths. Hence, an analogous 
procedure requires the ln s plane equivalent of a straight line, slope, 
and length. 

7 The ln s plane 

The s plane is defined by the rectangular coordinates s= Vt yw 
and the polar coordinates s = /s[ Ze . The ln s plane, hereinafter 
denoted by w, shall be defined by rectangular coordinates w = u+ jv. 
The conformal mapping relation between the planes is then 
w=u+ jv = In(s) = In /s/ + j@ . Hence, equating reals and imagin- 


aries: 
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Thus the ordinate of the w plane measures © and the abscissa 
measures In |s| . Referring to Fig. 1, the w plane, note that © extends 
from -00 to +00, corresponding to more than one polar coordinate revolu- 
tion in the s plane. Therefore, every point on the s plane is repre- 
sented, in one to one correspondence in the w plane, in the infinite 
strip bounded by O<6 ae 2 . Hence, this strip provides single 
valued correspondence. 

The scales of © and ln |s| are equal. This choice was made to 
retain the properties of conformal mapping. The abscissa is linearly 
calibrated in ln | s | with one unit of ln |s| equal in length to one 
radian of the ordinate. Then a logarithmic scale division calibrates the 
abscissa directly in /s/ . This is a most important point; for by using 
four-cycle semi-log paper, any four decades of © vs js} may be plotted. 
The cross-hatched area on Fig. 1 depicts such an area, and corresponds 
to the cross-hatched area on the s plane shown in Fig. 2. 

Fig. 3 shows this same area on a full-sized sheet of semi-log paper. 
The horizontal lines (© = constant) represent radial lines through the 
origin of the s plane; the vertical lines ( /s/ = constant) must repre- 
sent circles around the origin of the s plane. In the s plane, any radial 
line can be calibrated in length by the family of circles about the 
origin, j|s}| = 1, 2,...@, Hence, the horizontal lines in the w plane 
are calibrated in length by the vertical lines, |s/ = 1, 2,+42+ © 
Summarizing, the semi-log paper provides all the s plane radial lines, 
already calibrated in length. Further, the slope of these radial lines is 
the ordinate scale. 

To analyze the general line, not passing through the origin, examina- 


tion of the relation between Gartesian coordinates, T, Ww, and u, v, 
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shows that: 


w= In (S) 
WwW UtAV Uwe 
tise AbY ; 
ro po el a atlea engin 


Equating real and imaginary terms: 
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In (tT) = ut In. cos v In (ws) = UtIh = /aV 


These relations define families of transcendental curves represent- 
ing the vertical and horizontal lines of the s plane, Fig. 4 shows these 
families in the w plane, covering the same four decades as Figs. l, 2, 
and 3. The w plane vertical and horizontal lines have been omitted for 
clarity, but the scales are present. Selected members of the family have 
been labeled in "s plane counterpart" notation for easy reference, 

Notice that the form of all the curves is identical. Examine the 
three curves: J =-0.5, Vo= a VY =-2, Change in the value of 
VY shifts the curve laterally. Reference to the defining equation shows 
this effect. Here, ln W is an arbitrary constant. Changing its value 
changes the value of u, equally, throughout the length of the curve. 


Next examine the curves wW—-+/ and w=-+2Z. The same lateral 


shift occurs as between VY = —| and v= —Z.. This observation is 
expected from the defining equation Ir wm A F In sin Vv and 
hence, K= In w — In sin Y . This effect is generalized 
below. 


At this point it is well to note that a term ln Ywith ‘\~ taking on 





negative values is really not a cause for concern. In the s plane, 
In(+y—) is a positive distance measured along the ©& = 0 line and 
In(- v ) is a positive distance measured along the @ = a1180° line, pro- 
viding the correct minus sign. 

Before moving to the general line, Gore = ay = on 
some helpful geometric properties should be examined. In the s plane, 
the family of lines W=-@...-/ 04/2, ,.. $@ 
are all symmetric about the line Y= 0 which is the ae axis. The 
foailyjof Vines y= —COm. “ipo tl) 2 . a) lee are 
all symmetric about the line W = Q, the horizontal axis. Thus a line 
of symmetry is a line drawn from the origin, perpendicular to the line 
in question. Since angles are preserved in conformal mapping, we expect 
these same families of lines (e.g. VY = %-/,:.—%) to be symmetrical 
about their respective perpendicular radials (e.g. U/ = 0) in the w plane. 
Reference to Fig. 4 provides contactors 
3. The general line 

Let the general line (not through the origin) in the s plane be 
fess fr) b . Fig. 5 shows such a line. Define the 
perpendicular from the origin to this line: "the perpendicular radial," 
and the length of the perpendicular radial ''d,"" The radial parallel to 
W= ma+b (i.e. W='nT ) is defined: "the parallel 


radial." 


Inote: The line (W = 0) is the same line as (6 = 0, 'S0"7 360°) 
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but satay = Lan “A, where “ = angle of inclination of line 
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e (sin V COS A-~ SINK COS v/ — bcos x 
, = b cos @ 
e (sin(v-a) = 5 


LA - lr sin (Vv ry, pote |r (b cosa) 


But Ut + lh siny a lh ee is the transformation equation 
for the Cartesian coordinate lines. Hence: the transcendental curves 
corresponding to all the straight lines in the s plane, not passing 
through the origin, have the same form as the Cartesian coordinate 
(vertical and horizontal) lines. Hence, the shape is the same, and they 
differ only in positioning. Therefore, one template may be used to draw 
any non-radial straight line provided adequate rules for positioning it 
are developed. 

Fig. 6 illustrates such a template. The template represents every 
non-radial straight line, mapped onto w, just as surely as the edge of a 
triangle or T-square represents every straight line in the s plane. 

For clarity, rules for positioning the template, reading slope, and 
reading length are listed now, with proofs to follow in the next section: 


a. The template is always used in the same orientation, with the 
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line of symmetry horizontal. (Place template on a T-square and 
position vertically and horizontally, keeping the orientation the 
same, ) 

b. The slope of any line corresponds to the template slope lines, 
read on the ordinate scale of the w plane. There are two such slope 
lines, always 180 degrees apart. These correspond to “sense'' on the 
line. 

c. To draw any line, knowing one point and slope, slide template 
vertically until slope lines match required slope on ordinate; then, 
slide horizontally until known point is on the template. 

d. To connect any two points with a "straight line:" plot these 
points on the w plane Ge. /s / ). Then slide vertically and hori- 
zontally until both points are on the template, and draw the line. 
(Line of symmetry and slope lines must remain horizontal--you are 
positioning, not orienting the template.) The template will only 
coincide with the two points at one place, since there is only one 
straight line between two points. The slope of this line can be 
read at the slope lines. 

e. Where “sense” is important, the 180 degree slope ambiguity is 
solved as follows: When connecting point one to point two, follow 
the template line from point one to point two and continue in that 
direction to the correct slope line. This line then reads the slope 
from point one to point two. 

£. The template is calibrated for length. The distance between any 
two points on the line is equal to (P, - P.dK, where Py and P, are 


the template scale readings at point one and point two respectively: 





and K is a scale multiplier, varying with the length of the per- 

pendicular radial. Values of K can be read directly on the |s/] 

scale at the template arrow. 

Using this template, straight lines, slopes, and lengths can be 
plotted as quickly on the w plane as a ruler and protracter plot them 
on the s plane. 

4. Justification (casual reader may skip) 


The w plane equation of the general straight line is 
iu + ee CV = /n (bcos x) 


u-alnlb cos e In (s inva) 


Recall: 40 = /h Ww - |p (s1n vy) was an equation of cartesian co- 

ordinates, and changes in (J shifted the curve laterally only. In the 

general equation, i ( b cos) corresponds to ln W, (band 

are constants for a given line.) Changing b will have the same qualita- 

tive effect as changing UU , i.e. shifting the curve laterally. 
Referring to Fig. 5, note b cos A= ad. Hence: 

K= d— /p (s rr ae 4) It has just been shown how, by 

moving b, d is changed and the curve is shifted horizontally. Now, hold 

d fixed and change the slope (or & ). Recall the Cartesian defining 

equations for vertical and horizontal lines: 

we Ing —/h (cosy) and 

oe [hw — [hn (Sin v) = I Ww — /h (cos(v-)) 

When d is held fixed and these two are compared, Vo = WwW, and the 

only difference is the CG term A= OO Tana IA = a Yo 

And it has been previously shown in Fig. 4 that this corresponded to a 


vertical shift of 90 degrees only. Hence, varying any parameter shifts 
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the curves position only. 

Reference to Fig. 5 shows that as |s | increases to large values, 

the slope of s approaches the slope of the general line. Since we plot 

|s} vs © in w, we expect the general line to become asymptotic to 
the parallel radial at large values. The parallel radials are 90 degrees 
(looking both directions) on either side from the perpendicular radial. 
And since the slope of the general line is equal to the slope of the 
parallel radial, we need only to read the slope of the parallel radial. 
Since all radials are the horizontal lines of the w plane semi-log paper, 
the ordinate scale is used. 

Consider the concept of length in the s plane. Given the correct 
length calibration of all lines passing through the origin, it is possible 
to calibrate any other line. Draw the parallel radial W=/»VY corre- 
sponding to the general line w=amot bh. Then mark off equal lengths 
on the parallel radial (-@...-2 -|,4+/ 42, .., © ) and erect 
perpendiculars to the parallel radial at these equally spaced points. 
This family of perpendiculars will then calibrate w=-=my+ 5 = in 
length. Now note that the semi-log graph paper provides all the radials 
pre-calibrated for distance in the w plane. Then, just as the family of 
lines an ra Vv =. h » with b the generating parameter, 
will calibrate w= Ag+ Jb in the s plane; so will this family, 
mapped into the w plane, calibrate hae 5 » Mapped into 
the w plane. 

Since the spacing between the calibrating family is not constant, but 
changes with ln |s } , length is not as simple in the w plane as in the 
s plane. 


A single calibrated template is desired, not a calibrating family of 
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curves. Hence, the template may be arbitrarily calibrated in length 
corresponding to length along the line Y= ef , using the family 
W = -0O.,,~2, ~!, O, +1,42,,,,.4@ 5 and a correction factor 
(multiplier) is used for other lines. 

In both the s plane and the w plane, the general line is orthogonal 
to and symmetric about the perpendicular radial, and the calibrating 
family of lines is orthogonal to and symmetric about the parallel radial. 
This fact is invariant with changes in slope. Hence, for a given d (the 
length of the perpendicular radial) the general line and its calibrating 
family and their corresponding lines of symmetry are always 90 degrees 
apart. On the w plane, Fig. 4, it can be seen that change in slope, with 
fixed d, merely moves the general line and its calibrating family up and 
down. There is no relative movement between the two. Hence, for a given 
d, one calibration is valid for any slope. d= 1 has been chosen for 
calibration of the template in Fig. 6. 

Thus, to make the calibration valid for any non-radial line, only cor- 
rection for changes in d is needed. As parameter d changes 
Elm of = |r (Sin y= 0 ) ]; the general line moves laterally, 
while the calibrating family remains fixed--changing the calibration. 

The scaling function results from two important properties of the w 
plane; 

a. All radials are scaled logarithmically to read in distance. 

Hence, the horizontal distance from one to two is equal to the dis- 

tance from five to ten; or a horizontal distance d is a constant 

multiplier, regardless of position on the radial. 

b. The horizontal distance between any two members of the. family of 


calibrating curves is constant regardless of slope; i.e. dele | 
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and iJ—=2 have the same horizontal separation in the w plane at 
any V 
Refer to the defining equations; 


iL Se = (5 EN 


ee Ue = (ip ou dees v)- Un W_ —In sin ) 


A, ~U, = lh iw, ln a (not a function of v) 


Hence: as the general line moves horizontally with changing d, each 
point on the line crosses an equal number of calibrating curves. Hence, 
length calibration changes by an equal amount anywhere along the line. 

In the limit, the general line coincides with the parallel radial. 
Hence, horizontal movement (increase in d) by one octave will change the 
value of the calibration of the general line's extremities by one octave; 
hence, changes the entire line calibration by one octave. This movement 
corresponds to a change of d from one to two. Therefore, d can be used 
as a scale calibration multiplier for the whole line. dis defined as the 
distance from origin to vertex in the w plane, hence can be read at the 
arrow on the template. 
eve Plotting root loci 

Root loci may now be drawn in the w plane, Each sketching aid in the 
s plane is based on straight lines and their slopes and lengths. All 
these are now available in the w plane. 

Logarithmic root loci must be plotted on semi-log paper with a 
specific relation between ordinate and abscissa scale divisions, to pre- 
serve conformal mapping: namely, one cycle equals 131.8 degrees. Such 
semi-log paper is commercially available through Electro-Scientific 


Industries Corpn., 7524 Southwest Macadam Avenue, Portland 19, Oregon; 
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designated as ESIAC Form 44. This paper is fourccuale, two and one-half 
inches per cycle by 6.83 inches per 360 degrees, with five degree divi- 
sions and 30-degree heavy divisions. 
Plot root loci as follows: 
a. Choose and label the four abscissa ( | s| ) decades of interest, 
eng.4.,01, 0.1, 1,. 10,106. 
b. Label ordinates, zero to 360 degrees, 
c. Convert poles and zeros to magnitude and angle form, and plot. 
(Poles and zeros at origin can not, and need not be plotted.) 
d. Sketch root locus segments on real axis (6 = 0, 180°, 360°). 
(Start with pole or zero most positive and trace the real axis from 
that location toward minus infinity. Those segments of the real axis 
which are on the negative side of an odd number of poles (or zeros 
plus poles) are root loci.) Care must be taken in this matter since 
this means tracing left on the positive real axis and right on the 
negative real axis. Furthermore, poles at the origin must be 
counted when crossing from positive to negative real axis. When k 
is negative, substitute "even" for "odd" in this rule. 
e, Count poles, 1 , and count zeros, n. Then len root locus seg- 
ments extend to infinity. 
£f, Locate and plot centroid of asymptotes on the real axis, 
~- _ (Efes-Zzers) 
—— ae 
g. Angles of asymptotes are found as follows: 


rei) 
(1) For positive k: Angles = (2 . ne 
me / 2/ - / Ziv (2-nl-!) 


(2) For negative k: Angles --~- a (ian madians 


je-n] 


ll 





To plot asymptotes, place template slope line coincident with graph 
paper horizontal line, © = _ required slope, slide horizontally 
until template coincides with YW and draw asymptote, Repeat until 
all asymptotes are plotted. 
h. Angle of emergence from a complex pole may be calculated just as 
in linear theory. Keeping template slope lines horizontal, slide 
template until it connects the complex pole/zero with another pole/ 
zero. Read slope of connecting line at template slope line and 
record. Repeat for all poles and zeros. Then angle of emergence 
equals 180 degrees - S& (recorded slopes) for positive k, and 360 
degrees - S& (recorded slopes) for negative k. Place template slope 
line on this calculated value, slide over to complex pole/zero, and 
draw a short segment from pole or zero in template direction. This 
segment is the emergence from the pole. 
i. Sketch in locus. At any questionable points where accuracy is 
required, use the template as a Spirule; i.e., assume a trial point; 
measure angles from all poles and zeros (as is shown above) with the 
template. Then, if the sum of the angles equals 180 degrees (for 
positive k, or 360 degrees for negative k), the point is on the 
locus. If not, another trial is needed. Direction of new trial 
point is indicated by whether the sum of the angles is greater than 
or less than 180 degrees. 
j.- Root and gain point location on the locus is analogous to linear 
methods; i.e., k = [s+ of the loop transfer 

IT is +z 
function. Use template to measure distance from each pole/zero to 
point on locus in question and apply above equation. (Connect two 


points with template; length = (P, = PK where {= | s| as read at 
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II 
APPLICATION TO SENSITIVITY DESIGN 

l. Sensitivity design in the s plane 

Gain and root sensitivity are important considerations in control 
system design and analysis. An exact method with no approximations has 
been recently developed by Rung and Thaler (2). This method uses root 
locus plot in the s plane which is adaptable to the ln s plane. The bal- 
ance of this section is devoted to explanation of the method, defini- 
tions, and adaptation to the w plane. Reference (2) contains the proofs 
of the relationships between definitions to follow. 


Sensitivities are defined as vector quantities. 


——— 


| Lei GIS ah. yor S99 








Eq. 2.1 -- Root-Gain Sensitivity Y = Si /i- all i | 

} AS [SFR 
where q; is the root in question and L is loop gain. Thus Gain Sensi- 
tivity = change in root position per percent change in gain. (Since 
parameter changes may move the root in question in different directions, 


cme pees ae 


sensitivities must be vector a 





Eq. 2.2 -- Root-Zero Sensitivity 4 ae Z 7, ~ 
me a : ‘ qu 
Eq. 2.3 =- Root-Pole Sensitivity a — a se 


Pole and zero sensitivities = change in root ae per change in pole 
or zero position. Rung and Thaler (2) also show that the sum of root and 


pole sensitivities equals one. 


Eq. 2.4 -- A oy 7 z on, = = Using this fact along 


with equations 2.2 and >. vector representations of sensitivities are 


shown in Fig. 7. 
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Fig. 7 is a root locus plot and sensitivity vector diagram of the 
a Shae aie 
open loop transfer function Gis) —annigea neg pay 7i(s +3) 


in the s plane; k was chosen equal to 2.42, corresponding to a root at Fa 
as shown. 
The following construction procedure was used, Draw a vector from $4 
toward each pole with length equal to (1/distance from $4 to pole). 
ye 


i 4 ad 
These vectors define pole sensitivities, Sp, 9 Sp, 5 . 
2 


Draw a vector from Q, away from each zero, defining zero sensitivities, 
—_ _—y 


1 4 arr | | i | 
Sz, s Sze : Form the vector sum of Sp, ) Sp, 5 Sp, y Ser De, 


4 


_ CEE CLS ae ee 


ag a een ees 
and zero sensitivity vectors Po, Sp, at Sa, ,9 Phase is 


22. — 

—— / 
measured clockwise from U. For example, Cc a length of oP: ‘ x 

eae lage Scapa 
o Ens oO 
or Sp, =. 585276 J 
A _ 

Gain sensitivity Sh is directionally coincident with U which is always 
tangent to the root locus. (When gain changes, root must move along root 


é 


A. 
locus.) Magnitude of the S; vector is found from the relationship: 
Tae | aad, 








A S k a ze 
, ie or = Se, (Fu are 
“fT FR pag ee’ 


A — 7 
When the pole is chosen at the origin, =) ee Ye (+ Z.) 
eo 


k 


In this example —F, 2 | L228 Z. ge 








= | , o a ee 
or Sa ee Za S,=2 -—Liez 4 L326 
/ ) To Ze 
——— ae 








i (4 L.320°) (hez5£ 23°) —_ ny Ay rene ae 


k 


=) 
(measured counterclockwise as U is measured). 
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Since sensitivities are measured relative ue WU, the designer can 
manipulate the poles and zeros of compensation devices to achieve optimum 
placement of U, in addition to stability and response considerations. 

This technique is demonstrated clearly by Rung and Thaler (2). 
pe sensitivity design in the ln s plane 

For use in the w plane, this method is most practical and rapid using 
a separate graph for vector diagrams. Since straight line segments of 
given length and slope can be drawn with the template, vector diagrams 
can be drawn on the w plane, but the method is slower and gains no advant- 
age. After root locus and root position (or desired root position) q; 
have been plotted, use the template to measure and record distance and 
phase angle from q; to each open loop pole and zero. Then the vector 
diagram can be constructed on separate linear coordinates. The following 
example problem is taken from Rung and Thaler (2), exactly, but performed 
in the w plane to demonstrate the entire method. 

3. Design example 


The plant to be compensated has the open loop transfer function 


where 


Gis) = SCS 35t)) ~ S084 1,43) (St 7.39 
Kee diffe 


Dynamic bandwidth requirements lead to the desired location of the domin- 
ant roots at -q = -0.2¢ 50.35 = 0.4 120. that is, q = 0.4/-60. All 
three plant poles are subject to fluctuations, It is desired to design a 
geeetde compensator, satisfying the above dominant root requirement, and 
in addition, guaranteeing a minimum value of sensitivity of q to the poles' 
fluctuations. Fig. 8 shows plant singularities, R ) ei ’ Fe, and 


desired root location, qe Using the template, measure and record distance 
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and phase angle from each singularity to the desired root position. This 
information will first, determine the phase angle needed in the compen- 
sator to force the root locus through the desired root position, and 


secondly allow the vector diagram to be drawn. From Fig. 8: 


Singularity Phase Distance 1/Distance 
Po 120 0.4 PD 
Py 16 Teel 0.787 
P, = 56? 3.09 0.323 
142 


Phase angle needed from compensator is thus -180° + 142° = «38s indicat- 
ing a phase lag filter. Now we may draw the vector diagram (Fig. 9) 
assuming qs is on the root locus (remembering vectors are drawn from roots 
toward poles). Note that the ccupensa tom wanieess change the pole sensi- 
tivity vectors, but will rather add to their vector sum and change U, the 
reference, On Fig. 9 the vector sum of the sensitivity vectors is drawn, 
Or. Qt is the uncompensated U' vector. The sensitivity vectors result- 
ing from the pole and zero of the compensator will vectorially add vont". 
and generate the locus of U as they are moved along the -jW axis 

(always with correct separation to provide -38° phase lag). Rung and 
Thaler (2) show this locus to be a circle with center at I, the tip of 

the "uncompensated U vector," with radius r = — 217 d ; where () = 38° 
(the phase shift of the compensator), and WW = 0.35 (the imaginary part 
of the dominant root q;)- 

Stepwise: the compensator pole will have a sensitivity vector of length 
1/QP where QP is the distance from pole P to Q. As P moves along the 
negative real axis, 1/QP varies from a maximum when P is perpendicularly 


below Q in the s plane to a minimum when P is at infinity. Added 
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eS 


vectorially to in 1/QP generates a circle of radius R = 1/2W,= 1/0.7 = 1.43, 

with center vertically below I. This circle is drawn on Fig. 9 and labeled 

the 'M'’ circle. Then the sensitivity vector generated by the compensator 
—_— 

zero, 1/QZ, is added vectorially to the "M" circle, resulting in the U 

circle, with paw 3/4 2D, center at I (proof in reference two). The 

U circle is drawn and labeled in Fig. 9. 

In this example, minimum pole sensitivity has been specified, hence the 
longest possible U vector is desired. QU, the intersection of the U circle 
with the extension of QL, is this vector (desired U). Now, to find P and Z 
from U, proceed as follows; Draw OJ perpendicular to QU. Measure arcs 
JIN= M= ~ . Then, the direction of IN is the direction of QZ (i.e., 208°), 
and the direction of ™ is the direction of OP Gives, 246°). (Proof in ref- 
erence two.) 

Using these directions and the template, P and Z are located on the 
~180° axis in Fig. 8. (P = 0.9/-180° and Z = 0.36/-180..) 

4, Summary 

This method has designed the compensator to force the dominant roots to 

the desired point and ensured minimum root sensitivity for cascade compensa- 


tion. Furthermore, the sensitivity of each pole is known. Sensitivities 


are tabulated below. a 
f e 
a, = 27,2" | ;ge Far 
Sp Soc) 2 OMe 247 Cee Sone ) 
: ; a“ Of} sn fO% 
ee = 0./724775° B= o 093 LES* (Lompense to) 
a ) | 
Sy. elo SS x Late s* — 0,2¢94FL 273 = SPe gd 
a. 


l 


Ti, AW hae) AA 


CY ov 
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Sensitivity design is obviously simpler in the s plane than in the 

w plane, But when a root locus can not be handled in the s plane, due to 
widely spread open loop poles and zeros, or when desired root location is 
very near the origin, this w plane technique does work. There is no loss 
in usefulness or generality caused by plotting the vector diagram on 
linear coordinates. Choose the scale so that the largest sensitivity 
vectors plot easily. Then any vectors which are too small to plot 
accurately can be ignored, since their effect on the UV’ vector will be 


negligible, and small sensitivities are desirable. 
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III 
COMPENSATION DESIGN EXAMPLE 
Ie Compensation 
This example demonstrates: 
a. Clarity of w plane root locus analysis for compensator design 
where a pole/zero pair is very close to the origin and other poles 
and zeros are far removed. 
b. Relationship of frequency and transient response to w plane root 
locus. 
c. Breakaway~point calculation by plot of k versus real s (Wheeler 
Plot) directly on root locus plot. 
d, Graphical calculation of residues. 
e. Direct reading scale of damping factor of dominant roots on w 
plane root locus plot. 
A stabilization-type servo is designed as an example by Thaler and 
Brown (1). Compensation of the roll-stabilization loop in the w plane 


follows: 






R¢s) 






“SP 


FESS(ST SSS) 


The uncompensated open loop transfer function is: 


— _/35 09° aw 
G PPro S(st oss\St3azg 


A symmetrical compensator is chosen to meet the requirement of preserving 


the gain. Trial and error methods lead to a trial compensator of 
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(S+0.°5) (35 WD), 


G (Oofdls pad 7) 


The compensated system then has an open loop transfer function: 





G(s) = 


SCSt.055)(S+. 59S) S42). YCs +733) 


Fig. 10 is a root locus plot of the uncompensated system in the w 
plane, showing the closed loop complex roots corresponding to k = 135,000, 
and the instability of the system; the phase margin is zero. The locus 
breakaway point was determined by plotting gains versus real s. This plot 
can be made directly on the root locus plot, as shown. 

Fig. ll is the w plane plot of the root locus of the compensated 
System, plotted over eight decades. The superiority of the w plane is 
clearly demonstrated here, in its property of showing all the poles and 
zeros, regardless of separation, and the complete locus. 

The complex closed loop root corresponding to k = 135,000 is shown. 
The phase margin is 30 degrees, and the damping factor is 0.5. Since 
phase angle is a vertical coordinate in the w plane, the cosine may be 
plotted as a vertical scale, to read damping factor directly. 

2. Frequency response 

Fig. 12 and Fig. 13 show the Bode diagrams corresponding to Fig. 10 
and Fig. 11, respectively. Corner frequencies occur in the same position 
on the Bode diagram as the open loop poles and zeros appear in the root 
locus, The frequency response can be calculated directly from the root 
locus in the w plane, in a manner analogous to that in the s plane. How- 


ever, the standard approximate methods are much faster and should be used, 
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unless accurate detail is required. 
or Transient response 
The transient response may be calculated directly from the w plane 


root locus plot by calculation of the residues. 






elitr aa /3sseoo (st ss areas) 
RE) ~ [#Gls) ~ 1350S H sat SeAoTHst. Sohn! (St 333) 


Using roots found by Fig. 11 to factor the denominator, 


Ci) _ __pgsvoe( stesso ety) 
Kos) (St sas 42,42Z54334, Yotne ty loi e)(s + bm plore é] 


/ 
For a unit step input, K¢s) ees: and 


oe .) = lsso00(s+, ss)G4+2,/7) 
S (St. sere) Gt2. 42) St334. 2st Rb 1p Ve lS+ F.C ~ Vib) 


The inverse transform is: 
0 al -—Z, ae — ye 


Ete = Co +C,€é Lae + C3¢& rm 
nee (er, 
alice aan att als 4 oe a oe 


Evaluation of Co, Ce. is done on the root locus plot by evaluation 
of the residues, Wectors are drawn from all other roots and zeros to 
root s = <-2.42. Lengths and angles are recorded, and e, is evaluated 
from the equation below. (The root at -0.552 is effectively cancelled by 
the zero at -0.550 and could have been neglected.) This procedure is 


repeated for each root and all coefficients are thus evaluated. See below: 
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(When a root is on the real axis, the angular summation will equal I, 277, 


or a multiple thereof. In the above calculations, angles of zero or 2/7 


are omitted, and angles of 77 are indicated by a negative sign in front 


of the factor. Angular contributions from conjugate roots cancel, and 


conjugate factors were multipled together before entering above.) 


—~,sS26 2 G2 ; . 
} oe ~52h.22 
Cl) a. 7444+ 00386 € Pijzeese —.coxsze 
A257 C-76) 21etJe L257 CRC gh bLt 
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IV 
FACTORING A POLYNOMIAL 
L. 


One real root, one pair of complex-conjugate roots 


Any polynomial, regardless of order, may be factored by successive 
applications of the root locus plot. 


The following example demonstrates the factoring of a third-order case: 


Si4+ 529 4+ /0¢5 4260 


Set equal to zero and manipulate form: 
~ Coe e SP) em a 
S°-+ S52 5$*-+4+ (C45 4 2002C 
S* a S22 5S ~ + (eS = - zoe 


| » 4 
Sses25*¢/0985 





- 2 O00 


Se aie JE 


ee.  ~) ae 
5(5 Sfp eS ™" /o4) 





Ze 0 





“ ee, Zoe Seagir dae_ > ean 
S(st4sezs-10# 


S(SP4OAI(S bei1) 


Closed loop roots of this function are factors of the original polynomial, 
(see Fig. 14). Roots were found at 


— soLl¢ge’ = -sotgo 
- tao ee = —-/ Bs 
9 See 
by methods of Chapter I. Thus the factored form is: 


si4. S25 + /04S5 4200 (steols+4 west! gly 
2. Two pair of complex-conjugate roots 


Factor 


st 18S 74 7572-21084 25°0 
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Set equal to zero and manipulate, 

+ Ss 
Sgots + +75 32655 245250 
Stee es * +. G75" =m —-250 


= 2 = 21 Sa/1/¢ 
oe 7 ae ZA) (S+I~ 42,45) 


Closed loop roots of this function are the desired factors. 
To demonstrate the principle of successive root locus plots, an 
alternate manipulation will be used here in the solution of this fourth- 


order case. This technique can be applied to any order. 





g. . aaa 2. ae 
Se eS dee ee / OS Se 
- 2=°O bess, 2gC 
Bk. —/ = SPI S34 E785 “2/08 S(s*4/45*4 675-219) 


Now, to factor the cubic in the denominator, set it equal to zero and 


manipulate as before: S + 10S —F eyeeaee e 


—| Le Fe 
— SCs tHpiest 67) ~ SCS+ qt petystTge?? 


Closed loop roots of this function are determined by the root locus plot 


in Fig. 15. These roots are: 

§=2/7346° =/,73 + yo 

mee // 2 /53,5 = —e FI B42 

6 = NL 2b6,5°= -Y, tie TZ 

Thus the cubic is factored: (S—/,73/(S4+7 VIF ZF, 7H S47, righ 
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Then a root locus plot is made of Eq. A. 


-/ = ac 
SCS 1, 73)(St Tat C1 gto) ohne ge = 





Fig. 16 shows the closed loop roots to be: 

in 2. 7S = | ae 
e/a S a= of “Tage! 

S- ff! Lz20Gs =a =/Oe p- s7 


— 


s = // 2LIF3 5s —/O + ps 


ences of ae (FS ie 7 Gieao8 wees = (stor Ys roy VE-ltg )G-1-24 
For large-order polynomials, manipulations should be made to place 

the largest factorable (at least quadradic) group of terms in the numer- 

ator. Then, if the denominator is still unfactorable, the process is 

repeated. This procedure will reduce the polynomial size by three orders 


per plot. 
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V 
CONCLUSION AND APPENDICES 
ve Conclusion 

Root locus plots may be sketched or constructed in the ln s plane 
whenever the s plane plot fails to provide satisfactory resolution of the 
locus. Plotting in the ln s plane is somewhat more tedious and slower 
than plotting in the s plane, but does provide needed detail, while show- 
ing all the poles and zeros and all of the locus of interest. The 
s plane should be used in preference to the ln s plane whenever it is 
adequate. 

This paper lists a straightforward, step-by-step procedure for 
plotting in the ln s plane, and includes mathematical proof of the theory. 
The method is demonstrated in three practical examples in Chapters IT, 
III, and IV. 

Root locus plots in the ln s plane provide another useful tool in 


the design and analysis of feedback control systems. 
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